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We investigate rotational properties of a system of bosons with attractive interactions confined 
in a one-dimensional torus. Two kinds of ground states, uniform-density and bright-soliton states, 
are obtained analytically as functions of the strength of interaction and of the rotational frequency 
of the torus. The quantization of circulation appears in the uniform-density state, but disappears 
upon formation of the soliton. By comparison with the results of exact diagonalization of the many- 
body Hamiltonian, we show that the Bogoliubov theory is valid at absolute zero over a wide range 
of parameters. At finite temperature we employ the exact diagonalization method to examine how 
thermal fluctuations smear the plateaus of the quantized circulation. Finally, by rotating the system 
with an axisymmetry-breaking potential, we clarify the process in which the quantized circulation 
becomes thermodynamically stabilized. 
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I. INTRODUCTION 

Superfluidity represents the assembly of complex phe- 
nomena such as persistent currents, quantization of cir- 
culation, nonclassical rotational inertia, and topological 
excitations Q, 0- Gaseous Bose-Einstein condensates 
(BECs) offer a testing ground for these phenomena be- 
cause of their great flexibility to realize various experi- 
mental conditions. For instance, the Feshbach technique 
makes it possible to control the sign and strength of in- 
teractions. Furthermore, optical and magnetic traps of- 
fer ideal containers of BECs, in which microscopic surface 
rugosities that give rise to dissipation are either absent or 
can be manipulated as tunable parameters Q • These ex- 
periments have now become possible in low-dimensional 
systems 0,11111101 by tightening the confinement in 
one- or two- direction(s). Low-dimensional systems are 
simple theoretical models for st udy i ng vortices, persistent 
currents, and solitons H OH EE El El El El El El El 
For the case of attractive interactions, BECs do not col- 
lapse and form solitons in one dimension 0, H E| • In- 
terestingly, a bright soliton forms also in two dimen- 
sions if the strength of interaction is made to oscillate 
rapidly [T^ . 

In this paper we study a system of attractive bosons 
confined in a one-dimensional torus 19] under rotation. 
When the excitation in the radial direction is negligible, 
such a system is described by the Lieb-Liniger model [2£j 
with attractive interaction, in which a rotating term is 
added to the Hamiltonian and the periodic boundary 
condition in a finite system is explicitly taken into ac- 
count. When the rotational frequency of the container 
is increased, a Hartree-Fock theory [lj shows that the 
angular momentum of the system is found to exhibit 
plateaus of the quantized circulation [^J |51 HJ. We 
investigate the stability of quantized circulation by em- 
ploying the Gross-Pitaevskii mean- field theory (MFT), 
the Bogoliubov theory, and the exact dia gon alization of 
the many-body Hamiltonian. In Refs. 0E1, it is found 
that quantum fluctuations become significant near the 



boundary between the uniform and soliton phases in a 
non-rotating torus, and that the boundary is singular in 
the Bogoliubov approximation. We will show, however, 
that in the rotating torus quantum fluctuations are sig- 
nificant only in the immediate vicinity of a critical point 
at which a normalized rotational frequency of the con- 
tainer f2 is integral and the dimensionless strength of in- 
teraction 7 is equal to -1/2; for other values of f2 and 
7, there are no singularities in physical quantities at the 
phase boundary because the soliton can be formed with- 
out passing through the singular critical point. 

This paper is organized as follows. In Sec. [nj the 
ground-state wave functions of the system under rota- 
tion are derived analytically within the Gross-Pitaevskii 
MFT. In Sec. IIIII effects of quantum fluctuations on 
the quantized circulation are examined based on the Bo- 
goliubov theory and on the exact diagonalization of the 
many-body Hamiltonian. The results obtained by these 
two methods will be shown to agree very well, demon- 
strating the validity of the Bogoliubov theory. In Sec. II VI 
the circulation is calculated as a function of f2 at zero and 
finite temperatures based on the MFT and on the exact 
diagonalization method. In Sec. the response of the 
system, which is initially at rest, to a time-dependent 
axisymmetry-breaking potential is examined in order to 
clarify the process in which the quantized circulation be- 
comes thermodynamically stabilized. 



II. ANALYTIC SOLUTION OF THE 
ONE-DIMENSIONAL GROSS-PITAEVSKII 
EQUATION WITH A ROTATING DRIVE 

A. Hamiltonian for the system 

We consider a system of N identical bosons with mass 
M which are contained in a rotating torus of radius R and 
cross section S, where the angular frequency of rotation is 
2f2 in units of H/2MR 2 . The Hamiltonian for the system 
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in the rotating frame of reference is then given by 



p2tt 

t=\ d6 
Jo 



ft(0){L - Vl) 2 ^(9) + ^(9)^(9)^(9)^(9) 



(1) 



where L = —id/d9 is the angular-momentum operator, 
9 is the azimuthal angle, and U — &iraR/S characterizes 
the strength of interaction, where a is the s-wave scatter- 
ing length. Here and henceforth, the length, the energy, 
and the angular momentum are measured in units of R, 
h 2 /2MR 2 , and H, respectively. In Eq. (JTJ we include for 

convenience the kinetic energy 2 / c iddftip of the rigid 
body which is a constant and only shifts the zero of en- 
ergy- _ 

It is shown pill], by using the single- valuedness bound- 
ary condition of the many-body wave function that phys- 
ical quantities of the system described by the Hamilto- 
nian change periodically with respect to with the 
period of one. The phase is therefore characterized by a 
phase winding number 



J 



n + \ 



(2) 



where the symbol [x] expresses the maximum integer that 
does not exceed x, and by a continuous variable 



CO 



CI — J, 



(3) 



which is the angular frequency relative to </, and the 
range of u> is limited to —1/2 < uj < 1/2. In Fig. ^a), 
we show J and uj as functions of fl. 

We first seek the lowest-energy state of the one- 
dimensional Gross-Pitaevskii equation (GPE) in the ro- 
tating frame of reference, 



(L-Q) 2 + 2ir 1 \^(9)\ 2 tf,(6) = fi^(9), 



(4) 



where fi is the chemical potential, and a dimensionless 
parameter 

UN 
"27' 



7 ; 



(5) 



gives the ratio of the mean-field interaction energy to the 
zero-point kinetic energy. The condensate wave function 
ip(ff) is assumed to obey the periodic boundary condition 
ip(0) = V>(2tt), and is normalized as J " \ip(9)\ 2 d9 = 1. It 
is appropriate to assume the form of the solution as 



#) = V^e« (6) 
where p is the number density and up the phase. 

B. Uniform-density solution 

The stationary state which circulates on the ring with 
a uniform density is a plane wave, 



[i = uj +7. 



(7) 
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FIG. 1: (a) Phase winding number J and angular frequency 
(relative to J) uj = Q — J as functions of fi. (b) Ground- 
state phase diagram. The bold curve corresponds to the phase 
boundary determined by Eq. JTIJ: 7 - 2(fi - J) 2 + 1/2 = 0. 
The upper and lower regions of the bold curve correspond 
to the uniform-density phase and the bright-soliton phase, 
respectively. The winding number J takes on a constant inte- 
ger in between the vertical solid lines in Fig. Q^b) and changes 
by 1 when crossing them. On the solid lines, the ground-state 
wave function is described by the Jacobian cn-function which 
has one node. On the vertical dotted lines, the ground state 
wave function is described by the Jacobian dn-function. In 
between the horizontal dotted line and the bold curve, the 
uniform-density state is thermodynamically unstable, and be- 
low the horizontal dotted line, the uniform-density state is 
dynamically unstable. 



The stability of the ground state is determined by the 
sign of the lowest excitation energy. Diagonalizing the 
Bogoliubov-de Gennes (BdG) equations 



d 



'!- n 



[1 + 4tt7|V'| 5 



u„ + 2ir"/ijj 2 v n = X„u n 
v n + 27T7V>* 2 M n 



(8) 



we obtain the excitation energies A„ and the correspond- 
ing amplitudes u n , v n with positive norms as 



A n = y/n 2 (n 2 + 27) - 2nuj, 
u n = A/-+e i(J+ " )e , 



(9) 
(10) 
(11) 



where the normalization constants Af^ are determined 
from the orthonormality condition 



as 



[u n (0)u* m (0) - v n {6)v* m (6)}d8 = S n , m , (12) 



Af d 



\ 



1 

4tt 



n 2 + 7 



yV(n 2 + 27) 



± 1 



(13) 



For repulsive interactions 7 > 0, the ground state of 
Eq. J3J| always takes the form of the uniform-density so- 
lution (JJJ, since the lowest excitation energy A_i (for 
-1/2 < u < 0) or Ai (for < ui < 1/2) from it is posi- 
tive for all 7 and fl. However, for the attractive case, the 
first excitation energy Ai or A_i becomes zero at 



9 1 
7 - 2uj 2 + - = 0. 



(14) 



In Fig. ^b), we show the phase diagram of the ground 
states with respect to 7 and ft. The phase boundary 1(14^1 
of the ground states is represented by the bold curve. 
The uniform-density solution becomes thermodynami- 
cally unstable for —1/2 < 7 < 2ui 2 — 1/2 (in between 
the horizontal dotted line and the bold curve) because 
the first excitation energy becomes negative, while it be- 
comes dynamically unstable for 7 < —1/2 (in the lower 
region of the horizontal dotted line) because the first ex- 
citation energy acquires an imaginary part. The uniform- 
density state is thus stable only when 7 > 2lu 2 — 1/2 (i.e., 
in the upper region of the bold curve). 



C. Bright-soliton solution 

For 7 < 2uj 2 — 1/2 (in the lower region of the bold curve 
in Fig. [Jb)), the uniform-density state is either thermo- 
dynamically or dynamically unstable, and a soliton state 
becomes stable. In this subsection, we derive the ground- 
state soliton solution of Eq. J3J . The definitions of several 
kinds of elliptic integrals and elliptic functions (2|| used 
throughout this paper are summarized in appendix 1X1 to 
make the present paper self-contained. 

Substituting Eq. JBJ into the GPE |gj and taking the 
real and the imaginary parts, we obtain 



/i = ft 2 



/\2 



2$V + 2?T7p, 



^+V^-2fl^ =«>. 



Equation i|16|) is integrated to give 



(15) 
(16) 



Substituting this into Eq. I|15|) yields 

. ... ,r 2 

pp 



7r7p J + Vp - — ■ 



W 1 



(18) 



This equation can be rewritten in the form of an elliptic 
integral 



dO 



dp 



v /4 7r7 p3 _ 4 /i/9 2 + ±y p _ 4 W 2 



(19) 



which has formally the same solution as the case without 
rotating term p"?| and is given by 



p{9)=N 2 





)-6 ) 


ni^j — r\m,\ 









(20) 



where dn(u|m) is the Jacobian elliptic function, < r\ < 
1, and mi = 1 — m. A constant n is given below in 
Eq. (|23J) and the parameter m will be determined later 
in Eq. (|36(l . We denote the complete elliptic integrals 
of the first and the second kinds as K = K(m) and 
E = E(m), respectively. Since the soliton breaks the 
translational symmetry, the solution (|20|l contains an ar- 
bitrary parameter 9 which is regarded as the center of 
mass of the total particles. From the normalization con- 
dition, the normalization constant AT is determined as 
Af 2 = K/ [2it(E - rjmiK)]. Substituting p in Eq. (|20|l in 
Eq. HBJ, we obtain 



K 



n 2 f Q f'> 



= 1 



Af 2 = 

H = 

7] = 

W = 



where /, / c , /d are defined as 



K 1 



2n(E~f]m 1 K) ir 3 \-f\' 

1 
2V 2 

fd . fc 




f 
fc 

u 



2K 2 - 2KE 
2 mi K 2 - 2KE 
2KE + tt 2 7, 



7T 2 7, 



7T 2 7, 



(21) 

(22) 
(23) 

(24) 



(25) 
(26) 
(27) 



and sgn denotes the sign function such that sgn(w) = 
-1 when -1/2 < w < (i.e., J - 1/2 < fl < J), and 
sgn(w) = +1 when < u < 1/2, (i.e., J < ft < J + 1/2, 
see Fig.nja)). The cases of u> = 0,±l/2 will be treated 
separately below. 

From Eqs. I|17|) and l|2U|l . the phase part is given by 

<p(8) = fdS 



4 



which is integrated analytically using the relations (|A7|) 
and (|A10p to give 



ip{6) = S16- sgnH^ 1 !! ( n; — 



/ K9 




[ n; — 


rn^J 







(29) 



where II(n; u\m) is the elliptic integral of the third kind, 
and the constants n and 82 are given by 



1 — rjini 



2mK 2 



I 2/ 

Jc/d 



(30) 



Since the wave function is single- valued, the phase ip sat- 
isfies 



<p(2ir) - ip(0) = 2irJ. 



(31) 



This condition can be used to determine the value of m 
as follows. In the present case, the parameters m and 
n satisfy the relation m < n < 1, and then using the 
relations (|A14() - (|A17|I . the elliptic integral H(n;u\m) at 
9 = 2ir reduces to 



IL(n;2K\m) 



K+-6 2 {1-A (e\m)} , (32) 



arcsm < 



mi/' 



(33) 



where Ao is Heuman's Lambda function defined in terms 
of the incomplete elliptic integrals of the first F{e\m) and 
the second kinds E{e\m) as 



Anfelm) 



[KE(e\mi) - (K - E)F(e\m 1 )] . (34) 



The phase at 6 = 2tt thus becomes 
<p(2ir) = 2irSl — sgn(o;) 



'2/d/c 



/ 



+ 7T(l-Ao) 



(35) 



Using the notation u = SI — J instead of SI and J, con- 
dition H31JI then leads to 



2ttM = 



'2/d/c 



/ 



+ 7T(1-A ), 



(36) 



which determines the parameter m. This equation has a 
solution only when 7 < 2uj 2 — 1/2 and < \uj\ < 1/2, 
i.e., in the region delimited by adjacent vertical dotted 
and solid lines and by the bold curve in Fig. |T[a) . 

Next let us consider two special cases: to = (O is 
equal to integer J, i.e., on the vertical dotted lines in 
Fig. IHb)) and |u;| = 1/2 (SI is equal to half-integer J + 
1/2, i.e., on the vertical solid lines). The ground-state 
wave function given by Eqs. I|2(J|I and l|29|) is simplified 
when \lo\ = and \lu\ — 1/2 according to the limiting 
values of several parameters summarized in appendix iBl 



When lu — 0, the amplitude \fp(0) reduces to the Ja- 
cobian elliptic function dn(u|m), and the phase f(0) re- 
duces to J9, 



n 171 \ 7r 



/?? 



AJ9 



(37) 



Substituting this solution into the GPE Q yields the 
chemical potential 



K 2 

A 4 = 2~(1 + «ii)> 



and an equation 



f d = 2KB + 7T 2 7 - 0, 



(38) 



(39) 



which determines the parameter m only when 7 < —1/2. 
The dn-solution includes the ground state in the rest con- 
tainer with Q = 0. 

When \ui\ = 1/2, the amplitude is given by the Ja- 
cobian elliptic function cn(u|m), and <p(6) reduces to 
(J +1/2)0, 



m = 



Ik 2 ? 



7T 3 |7| 



cn 



K( 



(40) 



The chemical potential and an equation that determines 
the parameter m are obtained as 

/i = -^-(l-2mx), (41) 
f c = 2m 1 K 2 - 2KE - vr 2 7 = 0. (42) 

The equation 1)42(1 has a solution only when 7 < 0. The 
wave function (|40|) has a node at 9 = 9q+tt and the phase 
jumps by an amount of 7r at the node. With increasing 
Si adiabatically, the vortex enters the ring through this 
node. 



D. Ground-state properties 

To illustrate the SI dependence of the soliton solution, 
we present in Fig. [21 the densities and the phases for 7 = 
—0.55, where the integral constant of the phase is chosen 
so that (p(n) = 0. The density profiles of the solitons 
depend only on the relative angular frequency \u>\. As 
the strength of interaction is increased, the parameter m 
approaches unity (see Fig. llll in appendix lB|l for all uj. In 
the limit of m — > 1, both dn(u|m) and cn(u|ro) become 
sechu which is a soliton solution in infinite space (lt| . 

The ground-state energy £ per atom of the uniform- 
density state is 



7 

2' 



(43) 
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FIG. 2: Densities p (solid curves) and phases <p (dotted 
curves) for several values of angular frequency to and phase 
winding number J = fl — u) with 7 = —0.55. The density 
depends only on \lu\. The phase difference ip(9 + 2n) — f(8) is 
given by 2nJ in Figs, (a) and (b), and by 2ttJ + tt in Fig. (c) 
because at \lu\ = 0.5 the wave function has a node at which 
the phase jumps by an amount of n. 




FIG. 3: (a) Ground-state energy £ per atom, and (b) chem- 
ical potential fi. Both of them decrease monotonically with 
increasing I7I and are periodic with respect to Q. 



Equation (|44|l reduces to the energy of the dn-solution in 
the limit lo — > as 



£ = - 



K 2 [(1 + m x )E + m x K] 
3n 2 E ' 



(45) 



and to that of cn-solution in the limit |a->| — > 1/2 as 

K 2 [(1 - 2m{)E - mi(2 - 3 mi )K] 
~ 3TT 2 (E- mi K) ' 



(46) 



The ground-state energy per atom £ and the chemical 
potential \i are shown in Figs.|3Ja) and (b). For a given 
7, the ground-state energy £ reaches minima for integer 
VL and maxima for half- integer fi, and is smooth every- 
where. In the regime 7 < —1/2, the chemical potential 
\i becomes maximal for integer Q and minimal for half- 
integer Q, and has kinks at the phase boundaries given 
by 7-2(0- J) 2 + 1/2 = 0. The type of the phase transi- 
tion is the same as that of the non-rotating case j24[; at 
the phase boundary, (i) £ is smooth, the first derivative 
of £ with respect to 7 or has a kink, and the second 
derivative of £ has a jump, and (ii) /i has a kink, and the 
first derivative of /i has a jump. Detailed behaviors of 
these quantities near the phase boundary are discussed 
in appendix IO 



III. EFFECTS OF QUANTUM FLUCTUATIONS 

Without rotation, the quantum depletion diverges at 
7 = —1/2 and the first excitation energy obtained by the 
Bogoliubov theory becomes gapless at th at p oint, which 
requires a modification of the MFT [24l |25| . We inves- 
tigate here whether in the presence of rotation there is 
such a singular point at which the effects beyond the Bo- 
goliubov theory are significant. 

We evaluate the depletion of the condensate calculated 




and that of the soliton is 

K [3E - (1 + mi)K] 



7 . 



2K 2 [3E 2 - 2(1 + mi )KE + mi K 2 ] 

37T 4 7 



(44) 
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FIG. 4: Quantum depletion N' (Eq. (1491 ') at the phase bound- 
ary calculated by the Bogoliubov theory. 
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by the Bogoliubov theory as 



N' = 



E \M0)\ 2 de, 

n#0 



(47) 



where v n (9) is the hole amplitude in the BdG equa- 
tion JHJ. If N'/N is of the order of unity, the validity of 
the Bogoliubov theory is not ensured. Since the excita- 
tion in the uniform-density regime is contributed mainly 
from the excitation with quantum number 1, the deple- 
tion in the uniform-density regime becomes 



N' 



1+7 



1. (48) 
-1/2, this is rewritten 



N' 



1, 



(49) 



At the phase boundary 7 = 2ui 2 
as 

1 

which is shown in Fig. ^ The quantum depletion di- 
verges at the boundary u = 0,7 = —1/2. However, as 
\ui\ is increased, N' at the phase boundary decreases and 
the depletion becomes much less pronounced as shown 
in Fig. QJ This result is also inferred from Fig. ^b). 
The uniform-density state becomes dynamically unsta- 
ble below the horizontal dotted line 7 = —1/2, and this 
line touches the phase boundary (the bold curve) only at 
w = 0. 

We compare low-lying energy levels obtained by the 
Bogoliubov theory with those obtained by the exact di- 
agonalization of the many-body Hamiltonian QJ. In 
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FIG. 5: Bogoliubov spectra \ n for |a>| = (solid curve), 
|w| = 0.2 (dashed-and-dotted curves), and \uj\ — 0.5 (dot- 
ted curves). The (quasi) zero-energy levels in the soliton 
regime correspond to the Goldstone modes associated with 
the translation-symmetry breaking due to soliton formation. 
Open circles show the excitation energy £i(£o) — £o{£o) ob- 
tained by the exact diagonalization of the Hamiltonian for 
N — 500 with \ui\ = 0.2, which agrees with the Bogoliubov 
spectrum Ai + A_i (bold curve) even at the phase boundary 
| 7 |=0.42. 



Fig. [SJ we show the Bogoliubov spectra A„ obtained by 
the BdG equation (JHJ) as a function of the strength of 
interaction for the angular frequency \lu\ — 0,0.2, and 
0.5. The curves represent the excitation energies from 
the uniform-density state for | 7 | < — 2uj 2 + 1/2, and 
those from the soliton for | 7 | > — 2uj 2 + 1/2. All lev- 
els are continuous at the boundaries | 7 | = — 2ui 2 + 1/2, 
which indicates a smooth crossover between the uniform- 
density state and the soliton state. When u> — (solid 
curves), Ai and A_i are degenerate in the uniform-density 
regime, but they separate into two branches, the first 
excitation and a Goldstone mode, in the soliton regime 
(I7I > 0.5). The Goldstone mode, which is associated 
with the translation-symmetry breaking, boosts the soli- 
ton along the ring without increasing the energy. As \uj\ 
is increased from zero (dashed-and-dotted curves), the 
degeneracy is lifted by rotation. 

Next we calculate the low-lying energy levels by the 
exact diagonalization of the many-body Hamiltonian 
to see how they deviate from the Bogoliubov spectra 
near the phase boundary. The procedure of the diag- 
onalization is the same as that without rotation |24| . 
Denoting the number of atoms with angular momen- 
tum k as rifc, the plane- wave bases are prepared as 
\ni -i c ,...,m Q -i,ni ,ni 0+1 ,...,ni 0+ i c ) where Z ~ J is 
the angular momentum of the condensate, and l c is the 
cutoff. Within the subspace in which the particle number 
and the total angular momentum are conserved as 



i +i c 

E 

k—ln —l c 



n k = N, 



lo+lc 



(50) 



h — ln — / c 



we perform the diagonalization of the Hamiltonian 
which is rewritten as 



lC = ^{l-ti) 2 c\ci 



1 

2N 



E 

h i in 11 



m-\-n—k — l 



(51) 



We denote the ground-state angular momentum as Cq 
that gives the lowest energy eigenvalue, and the ground- 
state energy as £q(£o). Open circles in Fig.|S]show exci- 
tation energy £i(£o) — ^0(^0 ) obtained by the diagonal- 
ization of the Hamiltonian (|51|l for N = 500 and \uj\ = 0.2 
with cutoff l c = 1. The results agree very well with the 
Bogoliubov spectrum Ai + A_i (we represent here only 
the excitation that conserve the total angular momentum 
£0). We have also confirmed that the difference between 
the Bogoliubov spectrum and the exact one becomes even 
smaller as \u)\ is further away from zero, consistent with 
the analysis of the depletion of the condensate. The Bo- 
goliubov theory is thus vindicated except for 7 ~ —1/2 
and Q ~ J. 
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IV. QUANTIZED CIRCULATION AT ZERO 
AND FINITE TEMPERATURE 

A. Quantized circulation at zero temperature 

The magnetic flux is excluded from a superconductor 
when the applied magnetic field is below a critical value. 
The analog of this Meisner effect in superfluid system is 
the Hess-Fairbank effect |2l|. in which the system is not 
set into rotation when the frequency f2 of a rotating drive 
is below a critical value £l c . When O exceeds ft c , the cir- 
culation of the system, which is defined as the integral of 
the superfluid velocit y a l ong a closed contour, is quan- 
tized in units of h/M |22j,|23|, in analogy with the case of 
a type-II superconductor in which quantized vortices en- 
ter the system when the external magnetic field exceeds 
the lower critical field H Cl . The applied magnetic field 
and the magnetic flux in the superconductor correspond 
to the applied rotation and the angular momentum (L) 
of the superfluid, respectively. 

We calculate the angular momentum of the ground 
state (£> = N ^i^{9) {-ide)^{9)d9 1 where tjj is the 
mean- field solution JJJ) with Eqs. ifflTjl and obtained 
m Sec. EU and we find the analytical expression for the 
expectation value of L per atom as 

lt\ IN- / J ' 7 > 2a; 2 -1/2, , 

{L) /N- | j + U} + 2 nW, 7 <2 w 2 -l/2, (52) 

where W is given in Eq. (|24|l [see also Fig. HOf a)] . In the 
limits \w\ — * and |cj| — > 1/2, J + ui + 2irW reduces to J 
and J +1/2, respectively. Equation l|52() is shown in Fig.0 
as a function of £1 and 7, where the plateaus correspond 
to the uniform-density regime, and the crossover regions 
between plateaus correspond to the soliton regime. For 
7 < —1/2, the stable uniform-density state does not exist 
and the plateaus disappear. 

Since the quantum fluctuation is small, Eq. (|52|l cor- 
rectly describes the ground-state angular momentum ex- 
cept for 7 ~ —1/2 and f2 ~ J. Both Hartree-Fock the- 
ory 01 an d Monte Carlo calculation [27j show that this 
system exhibits the Hess-Fairbank effect in a certain pa- 
rameter regime, which is consistent with our results. 

We briefly comment on the same problem for repulsive 
interaction (7 > 0). According to the MFT, the ground 
state with repulsive interaction has a uniform density for 
all parameters. The expectation value of the angular mo- 
mentum then increases stepwise like the noninteracting 
case with 7 = in Fig. The repulsive bosons do not 
prefer the mixing of different angular-momentum states 
since it costs the Fock exchange energy. 

B. Quantized circulation at finite temperature 

We examine the effect of thermal fluctuations on the 
quantized circulation at finite temperature. The total 
angular momentum of the system is obtained from the 



(L)JN 




FIG. 6: Expectation value of the angular momentum per 
atom at zero temperature. The plateau regions correspond 
to 7 > 2(fi - J) 2 - 1/2. When 7 < -0.5, the plateaus disap- 
pear. 

derivative of free energy T with respect to Q as 

<^ = "~ + ^ (53) 

where 

T h 2 /(2MR 2 ) 1 1 ' 

is a dimensionless temperature with ks being the Boltz- 
mann constant, (■ • • ) T denotes the ensemble average at 
temperature r, and the second term in Eq. i|53[) corre- 
sponds to the angular momentum of the rigid body aris- 
ing from the constant term of the Hamiltonian Q . 

To evaluate the free energy, we employ the exact diag- 
onalization method of the many-body Hamiltonian. All 
low-lying levels £ n have been obtained in Sec. lIIII and we 
use them to calculate the free energy 

T= -rln^]e- £ " /T . (55) 

a 

Figure [7] shows the angular momentum for several tem- 
peratures calculated from Eqs. H53[) and (|55|l . where the 
mean-field result at r = is also presented for compari- 
son as dotted lines. At absolute zero, the result obtained 
by the exact diagonalization agrees well with the mean- 
field result, which supports the validity of the MFT con- 
sistently with the results in Sec. lIIII As the temperature 
is increased, thermal excitations wash out the edge of 
the step of the circulation, and the angular momentum 
of the system approaches that of the classical fluid, i.e., 
(L) T = NCI. 

In the present paper, we have focused on the prop- 
erties near r = in this calculation; however, in prin- 
ciple, this method can be extended to higher tempera- 
ture by increasing the cutoff angular momentum as long 
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0.02 



0.1 



0.2 0.3 

n 



0.4 0.5 



FIG. 7: Angular momentum per atom for 7 = —0.2 and 
N = 500 at temperatures r = 0,0.5,1, and 2 obtained by 
diagonalization of Hamiltonian l|51|l . Dotted lines show the 
result of the mean-field theory at r = 0. 



as the excitations of radial modes are negligible. Near 
|o;| ~ 1/2, though the diagonalization gives less accurate 
results than at uj ~ because a larger number of bases 
are needed in order to allow the state to have a node, 
the accuracy can also be improved by the increase of the 
cutoff angular momentum. 

Finally, we consider an experimental situation. To be 
concrete, let us consider 7 Li; then Eq. (|54J) leads to 



T~ 34 



{R [/mi]) 2 



[nK]. 



(56) 



Thus for a torus with R = 1 /im, r = 2 in Fig. [7| corre- 
sponds to T — 68 nK which can be achieved with current 
experimental techniques. A torus geometry may also be 
set up by the technique of microelectronic chips |28|. 



PREPARATION OF THE GROUND STATE 
BY A STIRRING POTENTIAL 



Rotation of the system can actually be driven by a po- 
tential that breaks the axisymmetry of the system. As a 
concrete example, we consider a time-dependent poten- 
tial 



V(6,t) =V coa(9-2nt), 



(57) 



which stirs the system with angular frequency 20, and 
we fix the amplitude of the potential as Vq = 1CP 3 . We 
take as an initial state if>(6, t = 0) = l/v2~7r with a fixed 
strength of interaction 7 = —1/4, and let the system 



(a) n=0.3 



p 0.01 7tfW^wV^-WVW^V\M^AMwMV^W^^ ^ 



0.35 



21)1) 



400 



600 



800 



1000 



g 0.02 




FIG. 8: Time evolutions of the angular momentum of the 
condensate (bold curves) and of the amplitude \ip(9o,t)\ at 
Oq = 7T (solid curves) driven by a time-dependent stirring 
potential £J3 with 7 = -1/4 for (a) Q = 0.3, (b) ft = 0.35, 
and (c) Q = Q C i- Note that the oscillations of the density 
and angular momentum are enhanced close to Q C r — v/l/8 ~ 
0.354. (d) Maximum amplitude of the wave function as a 
function of Q. 



evolve in time according to the GPE, 



de 2 



<KM)- (58) 



Figures |H1 (a)-(c) show the time evolutions of the an- 
gular momentum of the condensate and the amplitude 
\ip(0 = 7r, t) I . One of the phase boundaries for 7 = —1/4 
determined by Eq. (O is fl„ = ~ 0.354. As Q 

gets close to f2 cr , the stirring potential causes the dynam- 
ical instability, and a significant growth in density is seen 
in Fig.|Sfc). Figure IH^d) plots the maximum amplitude 
|^(#o)| as a mnc ti° n °f which shows the resonance 
induced near the phase boundary. However, the angu- 
lar momentum of the system oscillates without damping, 
and the system does not reach any stationary state. 

To achieve a stationary state, we must introduce en- 
ergy dissipation. We therefore study the time-evolution 
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FIG. 9: (a) Time evolution of the angular momentum of an 
initially uniform state with a damping constant V = 0.1 for 
7 = —1/4 and Vo = 10 -3 . (b) Expectation value of angular 
momentum per atom obtained by the MFT for 7 = —1/4. 



of the system according to a generalized GPE 29, 30], 



(i-T)^(6,t) = 



4 + Q 



V&cos0 + 2wj\^{6,t) | 5 



<KM),(59) 



where T is a phenomenological damping constant. Fig- 
ure (a) shows the time evolutions of the angular mo- 
mentum for several values of the stirring angular fre- 
quency with r = 0.1. The system acquires a finite an- 
gular momentum, and the magnitude of the absorbed 
angular momentum converges to that of the thermody- 
namically stable state as shown in Fig. EJb) . We note 
that as long as fl lies in the same plateau region, the 
angular momentum converges to the same integral value, 
demonstrating that the ground state derived in Sec. ^ 
can indeed be prepared by the time-dependent stirring 
potential in the presence of energy dissipation and that 
the circulations is indeed quantized in the thermodynam- 
ically stable state. 

From these results, it is concluded that the system 
reaches a thermodynamically stable state after both den- 
sity fluctuations and energy dissipation. The oscillations 
in the density is drastically enhanced near the resonant 
frequency (dynamical instability). However, the system 
can acquire a net angular momentum due to energy dis- 
sipation (Landau instability). Similar mechanisms were 
found in the vortex lattice formation Hal. 



VI. CONCLUSIONS 

In conclusion, we investigated the rotational proper- 
ties of one-dimensional bosons with attractive interac- 
tions confined in a rotating torus trap. 

We derived the ground state wave function analytically 
within the MFT as a function of the strength of interac- 
tion 7 and of the rotational frequency of the torus O. 



A uniform-density solution and a bright-soliton one are 
found to smoothly cross over to each other. The density 
of the soliton depends on the rotational frequency ui rel- 
ative to J and has a node atQ = J+ l/2 when a vortex 
enters the ring. 

In order to investigate the validity of the Bogoliubov 
theory, we compared the excitation spectrum obtained 
by the Bogoliubov theory with that obtained by the ex- 
act diagonalization of the many-body Hamiltonian, and 
found that the MFT well describes the ground state ex- 
cept for f2 ~ J, 7 ~ —1/2, where the phase boundary and 
the borderline of the onset of the dynamical instability 
coincide, and hence quantum fluctuations are significant. 

The angular momentum of the ground state is quan- 
tized with respect to il in the uniform-density regime, but 
it is not in the soliton-regime. The circulations at finite 
temperature are examined by the exact diagonalization 
method. 

To understand the process in which the system is set 
into rotation, we considered a time-dependent stirring 
potential which breaks the axisymmetry of the torus. 
The potential is shown to induce the dynamical instabil- 
ity, causing the density to oscillate resonantly and trig- 
gering the acquisition of the angular momentum by the 
system. With energy dissipation, the quantized circula- 
tion is found to become thermodynamically stabilized via 
the Landau instability. 
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APPENDIX A: ELLIPTIC INTEGRALS AND 
ELLIPTIC FUNCTIONS 

We summarize definitions and interrelations of some 
elliptic integrals and elliptic functions |2|j that are used 
in the present paper. 

Elliptic integrals of the first kind F, the second kind 
E, and the third kind Id are defined as 



F(4>\a) = 



cW 



y 1 - ••=• ■ ■ 



sin a sin 



E(4\a) = / dO\fi 



lT(n; 4>\a) 



sin 2 a sin 2 1 



d6 



(1 — n sin 2 9) sin 2 a sin 2 



(Al) 
(A2) 

:(A3) 



Using a parameter m = sin a, Jacobian elliptic functions 
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are defined by 



0.5 



sn(u|m) = sm<p, 
cn(u\m) = cos</>, 



(A4) 

(A5) 

dn(u|m) = *J 1 — m sin 2 (f>, (A6) 

and these functions are interrelated as 

dn 2 (u|m) — mi = mcn 2 (it|m) = m [l — sn 2 (w|m)] .(A7) 

Using the parameter set {u,m} instead of {</>, a}, the 
elliptic integrals are also expressed as 



F(u\m) = / dv = u, 
Jo 



E(u\m) — I dn (v\m)dv, 
Jo 

H(n;u\m) 



1 — nsn 2 (u|m) 



(A8) 
(A9) 
(A10) 



Elliptic integrals are said to be complete when <j> = 
7r/2, and are usually denoted as 



E(m) =E[ - m) . (All) 



K(m) = F ( - 



Complete elliptic integrals of the first and second kinds 
are expanded for \m\ < 1 as infinite series 



K(m) = - 



1+1 - I ro+| — ) m 



1-3 



1- 3-5 

2- 4-6 



E(m)-- 



2/ m I 2-4/ 3 I 2-4-6/ 5 



(A12) 



(A13) 



The elliptic integral of the third kind U(n; <f>\a) has other 
expressions depending on the relation between m and n. 
When m < n < 1 and <fi = tt/2, it reduces to 



n(n;|\a) = if (a) + |<5 2 {1 - A ( £ \a)} , (A14) 



(1 — n)(n — sin a) 



£ = arcsm 



1-n 



(A15) 
(A16) 



where A is Heuman's lambda function defined as 



A o (0\a) = - [K(a)E((t)\90 o - a) 



{K(a) - E{a)} F(0\9O° - a)} . (A17) 
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1 


n 


m 


1 


W 2 








5 2 


oo 


oo 


e 


?r/2 






l + 2 7 


VI + 2 7 /4tt 

^1/(1 + 27) 
tt/2 



TABLE I: Limiting values of various parameters of the soliton 
solution. The limit \uo\ — > is equivalent to the limit /d n — > 
or 7 — > —2KE/-K 2 — 0, and |w| — > 0.5 is equivalent to / cn — > 
or 7 — » — 2(KE— miK 2 ) / -k 2 — 0. The limit m — > corresponds 
to the uniform-density limit \^\ — > — 2oj 2 + 1/2 + 0. 



APPENDIX B: LIMITING BEHAVIORS OF THE 
SOLITON SOLUTIONS 



We consider the limits — >■ and |w| 
soliton solution for < \u)\ < 1/2 given by 



1/2 of the 



P{0) 



K 2 

7T 3 | 7 | 



dn 



K 



(0 



m l — 77TO1 



ip{9) = ne - sgn{uj)52 n ( n 







n; 




V ^ 





1 

2tt 



(U-fc-f), 



where 



/ = 2K 2 -2KE -tt 2 j 

fc 

k 



2 mi K 2 - 2KE - tt 2 7, 
2KE + 7T 2 7 . 



(Bl) 
(B2) 
(B3) 



(B4) 
(B5) 
(B6) 



The limiting values and behaviors of several parame- 
ters are summarized in Tableland Fig. 1101 respectively. 
In the limit \u\ — » (/d — > 0), the equations (|B1|) and 
(|B2|) continuously become those of the dn-solution since 
rj -> and W -> from Fig. [TOI and Table |l| In the 
limit \uj\ 1/2 (f c — > 0), where the parameters behave 
as 77 — * 1, TU — > 0, Eqs. (|B1I) and 1B2I) reproduce the 
cn-solution. It can easily be verified by setting /a = 
or f c = in Eq. I|B3|I that in these limits the chemical 
potential for < \u\ < 1/2 continuously approaches 



K 



1 



/* = -fj (1 - 2mi) = ^(/ d - /), M = 1/2.(B8) 

The first derivative of the ground-state energy with 
respect to has a kink at the phase boundary, which is 
verified by the relations (|52(l and l|53|l as, 



d£_ 



2 (n - (L) /N 

0, H = 0,1/2, 

-4irW, 0<|w|<l/2. 



(B9) 
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FIG. 10: (a) (b) (c) Behaviors of several parameters appeared 
in the soliton solution for 7 = —0.4 (dashed-and-dotted), 
7 = -0.5 (dashed), 7 = -0.6 (solid), and 7 = -0.7 (dot- 
ted curves), (d) Incomplete elliptic integral of the third kind 
II(n; K0/n\m) for |a>| = 0.4 as a function of 6. 



Figure ^2 shows the solution to calculated numerically 
by solving 



/d = 0, 



'2/d/c 



27TM =W ^^ + vr(l - A ), < \u\ < 1/2,(B10) 

/c = 0, 



H = o, 
x M < 
H = 1/2, 



where the line of intersection between the curves and the 
7-w plane corresponds to the phase boundary 7 — 2uj 2 + 
1/2 = 0. 



APPENDIX C: GROUND-STATE PROPERTIES 
NEAR THE PHASE BOUNDARY 

We investigate the continuity of the ground-state prop- 
erties at the phase boundary. From Fig. ^2 we see that 
the parameter m becomes zero for all soliton solutions in 
the limit I7I — > — 2w 2 + l/2 + 0. The uniform-density limit 
of the soliton solution is hence mathematically obtained 
by taking the limit of to — > 0. Let 6 be a positive small 
deviation of 7 from the value at the phase boundary for a 
fixed u>. For \oj\ = 0, 1/2, the parameter m, K(m), E(m), 
and hence the physical quantities near the phase bound- 
ary are expanded in terms of 5 as follows. 




FIG. 11: The solution m calculated numerically 
from Eq. ifMlOl . The dotted curve on the 01-7 plane 
represents the phase boundary 7 — 2lq 2 + 1/2 = 0, and 
there exists a unique solution m in the soliton regime 

It I > -2uj 2 + 1/2. 



For the dn-solution (u) = 0), the phase boundary is 
7 = —1/2 and let 7 = —1/2 — 5. Using the equation 
/d = 0, we obtain 

to = 8<5 1/2 - 32(5 + 89(5 3/2 - 20(W 2 + 0{5 5/2 ). (CI) 

The expansion formulae (|A12(1 and l|A13|) become 

K=^(l + 25 1 ' 2 + 5 + V/ 2 + U 2 ^j + 0(S 5 / 2 ), (C2) 

E = —(l — 25 1 ' 2 + 5,5 - -«5 3 / 2 + —S 2 ) + 0{5 5 ' 2 ). 
2 V 4 2 / 

(C3) 

Using these expansions, the ground-state energy and the 
chemical potential are expressed as functions of 5 instead 
of m near the phase boundary 7 = —1/2 — <5. The jump 
of the second derivative of the ground-state energy and 
the first derivative of the chemical potential are obtained 
as 

C - E'l = -4, Mdn - fi' u = -2, (C4) 

at 7 = —1/2 and uj = 0. 

For the cn-solution (\u>\ = 1/2) at 7 = —S, the expan- 
sions are obtained in a similar way as 



to = AS -6<5 2 + 0(S 3 ), 



(C5) 



K = l(l + 6+h 2 ) .Our'). (C(,i 
E = l(i- S +l5 2 )+0(S*). (07) 

The ground-state energy and the chemical potential are 
expanded in terms of S as 



/'c 



I _ M - — 



0(S 3 ) 



(08) 
(C9) 
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which are continuously connected with those of the and have discontinuous jumps at 7 = and \ui\ = 1/2 by 
uniform-density solution at 7 = and \u>\ = 1/2. How- the following amounts: 
ever, the first derivative of the ground-state energy and 
that of the chemical potential are given by 

C =-\-{ + 0(5% (CIO) C " K = ~4' M - ~ ^ = ~ 2- (C12) 

/4 n = -\-^ + 0(5% (Cll) 
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